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Abstract 

We present the first measurement of the atomic mass dependence of central 
H~ and H + production. It is measured using a sample of 22,459 H~'s and H + 's 
produced in collisions between a 250 GeV tt~ beam and targets of beryllium, 
aluminum, copper, and tungsten. The relative cross sections are fit to the 

1 



two parameter function aoA a , where A is the atomic mass. We measure 

a = 0.924 ± 0.020 ± 0.025, for Feynman-x in the range -0.09 < x F < 0.15. 
13.85.Ni, 12.38.Qk, 25.80.Ls 
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The atomic mass dependence of strong interaction cross sections with nuclear targets is 
sensitive to the behavior of hadrons and quarks inside nuclear matter. In addition, knowl- 
edge of this dependence is needed to compare cross section results from experiments using 
different target materials. Many atomic mass dependence measurements have been made 
||TJ. Nevertheless, little exists in the literature for central hyperon production. We report 
here, the first measurement of the atomic mass dependence of central S~ and S + production. 

The atomic mass dependence of cross sections is frequently parameterized as 

a(A) = a A a , (1) 

where A is the atomic mass of the target. By using four different target materials, we are 
able to check the applicability of this parameterization, as well as making a measurement of 
a:. 

A model of the nucleus as a totally absorbing sphere gives a value for a of 0.67. For 
absorption cross sections, a is a little higher. For example, Carroll et al. measured 
a = 0.755 ± 0.010 for a 280 GeV ir~ beam. If a cross section were simply proportional to 
the number of nucleons in the nucleus, a would be 1.00. Earlier, we reported a measured 
value for a of 1.00 ± 0.05 ± 0.02 for D meson production [^J. 

The apparatus in Fermilab experiment E769 has been previously described (see Q and 
references therein). The targets were 26 foils of Be, Al, Cu, and W with a total nuclear 
interaction length of 2%. The foils were simultaneously exposed to the 250 GeV beam to 
minimize errors associated with flux measurement. A differential Cerenkov counter was 
used to measure the beam content and reduce contamination from kaons. The final sample 
consisted of 95% %~ with a contamination of 3% K~ and 2% p. 

The elements of the spectrometer relevant to this analysis are 11 silicon microstrip planes 
(1-30 cm downstream of the targets), 35 drift chamber planes (150-1750 cm), 2 multiwire 
proportional chambers (130 cm, 180 cm), and 2 magnets (290 cm, 620 cm) for momentum 
measurement. The electromagnetic and hadronic calorimeters were used only for on-line 
event selection. The two threshold Cerenkov counters downstream of the target were not 
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used in this analysis. 

Tracks of charged particles were reconstructed using hits in the detector planes. The S^'s 
were reconstructed using only the three final tracks produced in the decays S~ — > A + and 
A — > 7T _ +p (charge conjugates are implied in this paragraph and the following paragraph). 
The analysis focuses on events where both decays occurred between the silicon microstrip 
planes and the drift chamber planes. The tracks used to reconstruct S~'s were required to 
have hits only in the drift chambers and proportional chambers, no hits being allowed in 
the silicon microstrip detectors. 

We applied further criteria to select S~'s from the data. The reconstructed A mass was 
required to be within 5.25 MeV of the known mass. The shortest distance between the two 
tracks used to reconstruct the neutral A was required to be less than 0.7 cm. The shortest 
distance between the A track and the other charged pion track was required to be less than 
0.66 cm. The angle between the S~ trajectory and the direction from the primary vertex to 
the S~ decay vertex was required to be less than 0.012 radians. There were requirements on 
the geometric locations of the three vertices and on the charges of the three decay tracks. 
These were the most significant selection criteria. 

Figure [TJ shows the invariant mass distribution for candidate S~'s and S + 's after the 
selection criteria were applied. The figure shows a strong signal over a linear background. 
The signal was determined using sideband subtraction, not by fitting a function to the mass 
distribution. This eliminates errors associated with the determination of the shape of the 
signal peak. 

The combined S~ and S + data signals, before acceptance corrections and weighting are 
applied, are as follows: 1980 ±55 from a minimum bias trigger (no requirement on transverse 
energy in the calorimeters), and 20479±187 from a trigger that required greater than roughly 
5.5 GeV of transverse energy in the calorimeters. The data acquisition rates for these triggers 
were controlled using prescalers which were set to record a specific fraction of events that 
passed trigger requirements. The signals from the two triggers are combined using weights 
based on the known prescaler settings. Typically, the minimum bias events with less than 
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5.5 GeV of tranverse energy have weights roughly 20 times larger than the events from the 
transverse energy trigger. This causes the statistical error on a to be dominated by the 
statistical errors on the events from the minimum bias trigger. 

A full detector simulation was used to calculate the acceptances for each material in 
narrow bins of xp (width 0.015). The acceptance calculation was repeated for each data 
set listed in Table [l[ Acceptances ranged from 3% to 12%. The simulation modeled the 
geometry of the detector, the primary interaction, secondary interactions, pair production, 
multiple scattering, detector plane efficiencies, and all analysis cuts. A total of 1.8 million 
simulated events were generated for the acceptance calculation. The statistical errors on the 
acceptances were much smaller than the statistical errors from the data. 

The dominant systematic error is due to the uncertainty in simulating the average number 
of charged particles per event (the multiplicity). The average multiplicity increases with the 
atomic mass of the target. Since the targets were arranged along the beam direction in 
order of decreasing atomic mass, events produced in the higher mass targets also suffered 
more pair production and secondary interactions than those produced in low mass targets. 
The overall effect is to reduce the acceptance for the higher mass targets, because the track 
reconstruction efficiency is lower in high multiplicity events. This systematic error was 
studied by varying the multiplicity of generated events in the simulation and by making 
comparisons between the data and the simulation. The systematic error on a related to 
multiplicity was estimated to be ±0.023. The systematic error was somewhat higher in 
events where the S~'s and S + 's had low transverse momenta (pt), because the track densities 
were higher nearer the beam. 

We studied several other sources of systematic error. Systematic errors in a associated 
with measurement of the thickness of the target foils, and errors associated with the location 
of the reconstructed primary vertex were each estimated to be roughly ±0.007. Other 
systematic errors were estimated to be even smaller. These include errors related to inelastic 
collisions in the target that attenuate the beam flux, simulated S~ and S + momentum 
distributions, detector geometry, signal determination, beam contamination, and trigger 
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biases. 

The values of a were determined by fitting the two parameter function OqA 01 to four 
data points: the relative Be, Al, Cu, and W cross sections. The fit parameter cro simply 
normalizes the function. Figure |2] shows one of the fits. The x 2 °f this fit is 1.35 with two 
degrees of freedom. The function fits our data well. 

Table H shows measured values of a for several data sets. The largest data set covers the 
region in xp where the acceptance is large enough to yield enough statistics for a meaningful 
result. This large data set is subdivided in several different ways. When calculated separately 
for S + and 5~, a is the same within errors. As a function of both pp and Xp, the dependence 
of a is consistent with being flat within errors, but there is a rise near xp = 0.00 and in 
the highest bin of pp. Note that the systematic errors are strongly correlated in different 
subsets of data. In evaluating trends as a function of Xp or pp, the statistical errors are 
more important than the systematic errors. 

It is useful to compare our results to other experimental measurements. In PH, the 
atomic mass dependence for production of 7r ± , K ± , p and p was reported. That experiment 
used a proton beam to study central production as a function of pp and beam energy. For 
example, at pp = 0.77 GeV with a 400 GeV beam, a was measured to be 0.91 ± 0.01 for 
production of 7r~ and 0.98 ±0.02 for production of both K~ and p. Our result for the atomic 
mass dependence of central H~ and S + production is similar in that a is also a little less 
than 1.00. In Ref. 0, the S° atomic mass dependence was reported for 0.2 < xp < 0.8 using 
a 400 GeV proton beam. The S° is the particle most similar to the H~. Figures 31 and 33 of 
Ref. show a as a function of pp and xf. These figures show the long established trends 
of a increasing as xp decreases towards 0.0 and a increasing as pp increases. The xp figure 
shows a increasing from less than 0.5 towards a value a little less than 1.00 as xp decreases 
from 0.8 to 0.0 (no data for xp < 0.2 in This is consistent with our measured value for 
a. 

In summary, we present the first measurement of the atomic mass dependence of the 
cross section for central S~ and H + production. We measure a = 0.924 ± 0.020 ± 0.025, 
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having found that A a is a good parameterization of the atomic mass dependence of the 
production cross section. Our data are consistent with the value of a approaching 1.0 as xp 
approaches 0.0. 

We gratefully acknowledge funding from the U.S. Department of Energy, the U.S. Na- 
tional Science Foundation, the U.S. National Science Foundation Graduate Fellowship Pro- 
gram, the Wayne State University High Energy Physics Initiative, the Brazilian Conselho 
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of Canada. 
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FIGURES 

FIG. 1. Invariant mass distribution for H~ and E + candidates. Mass is calculated using the 
known A mass, the known ir~ mass, and the measured momenta of the three tracks (assumed to 
be pir~7r~ or ptt + tt + ). Only candidates that pass the analysis cuts are included. Only statistical 
errors are shown. Events are weighted based on trigger prescalers. The average weight is 1.56. 

FIG. 2. Relative cross sections for E + and H~ production (xf between —0.09 and 0.15, summed 
over pr) as a function of the atomic mass of the target. The solid curve shows the fit of the function 
a = aoA a , where uo and a are the fit parameters. Only statistical errors are shown and used in 
the fit. 
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TABLES 



TABLE I. a with statistical then systematic errors. The overall data set contains E +, s and 
H^'s with xf between —0.09 and 0.15, summed over pr- The other lines show results from subsets 
of the overall data set with the additional selection specified in the first column. 


Data Set 


a in a = aoA a 


Overall 


0.924 ± 0.020 ± 0.025 


~ + only 
E~ only 


0.905 ± 0.028 ± 0.025 
0.939 ± 0.026 ± 0.025 


p T 0.0 to 0.5 GeV 
p T 0.5 to 1.0 GeV 
p T 1.0 to 1.5 GeV 


0.906 ± 0.041 ± 0.035 
0.913 ±0.028 ±0.022 
0.988 ± 0.041 ± 0.020 


x F -0.09 to -0.03 
x F -0.03 to 0.03 
x F 0.03 to 0.09 
x F 0.09 to 0.15 


0.881 ± 0.042 ± 0.025 
0.981 ± 0.029 ± 0.025 
0.910 ±0.034 ±0.025 
0.918 ± 0.056 ± 0.025 



10 



} t 



++ 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



1.29 1.3 1.31 1.32 1.33 1.34 1.35 



Mass (GeV) 



